Long-range order versus random-singlet phases in antiferromagnetic 
systems with quenched disorder 
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The stability of antiferromagnetic long-range order against quenched disorder is considered. A 
simple model of an antiferromagnet with a spatially varying Neel temperature is shown to possess 
a nontrivial fixed point corresponding to long-range order that is stable unless either the order 
parameter or the spatial dimensionality exceeds a critical value. The instability of this fixed point 
corresponds to the system entering a random-singlet phase. The stabilization of long-range order is 
due to quantum fluctuations, whose role in determining the phase diagram is discussed. 



a 

s 

I 

o 

> 

\o 

O 

ON 

O 
vn 

On 

a 

O 

o 



X 



Quantum antiferromagnetism (AFM) has experienced 
a surge of interest in recent years, both in efforts to ex- 
plain the magnetic properties of doped semiconductors 
1 1 1 , and in connection with high-Tc superconductivity 
1 2,^. In the former context in particular, the interplay 
between AFM and strong quenched disorder is an impor- 
tant issue. Bhatt and Lee have modeled the weakly 
doped, insulating regime of these systems at zero temper- 
ature (T = 0) by an ensemble of randomly distributed, 
AFM coupled Heisenberg spins with a very broad distri- 
bution of coupling constants J, and employed a numerical 
renormalization procedure ||]. While, with decreasing 
temperature, an increasing number of spin pairs freeze 
into inert singlets, they concluded that the remaining 
spins give essentially a free spin contribution to the mag- 
netic susceptibility. The net result is a 'random-singlet' 
(RS) phase, with a sub-Curie power-law behavior of the 
magnetic susceptibility. The quantum nature of the spins 
thus prevents the classically expected long-range order 
(LRO) of either AFM or spin glass type, and explains 
the experimentally observed absence of LRO. 

Bhatt and Fisher have applied similar ideas to the 
highly doped metallic regime. These authors argue that 
rare fluctuations in the random potential always provide 
traps for single electrons, which then act as randomly 
distributed local moments (LM) to which the method of 
Bhatt and Lee can be applied, but now in a metallic 
environment. Their conclusion was that the LM cannot 
be quenched by either the Kondo effect or the conduction 
electron induced RKKY interaction. This leads again to 
a RS phase with a magnetic susceptibility that diverges 
as T — > 0, albeit slower than any power. 

These results raise the important question whether, 
and how, AFM LRO can ever exist in a disordered sys- 
tem. Intuitively, one expects quantum fluctuations to 
weaken the metallic RS phase since they enhance the in- 
teraction of the isolated electrons with their environment. 



One should therefore wonder whether quantum fluctua- 
tions can restore LRO by suppressing the RS phase that 
would otherwise pre-empt an AFM transition. 

In this Letter we address these questions by studying a 
model of an itinerant AFM with a spatially random Neel 
temperature. We first describe our main results. We find 
that quantum fluctuations do indeed restore a LRO AFM 
phase, provided that the order parameter dimensionality, 
p, is smaller than a critical value Pc which depends on the 
spatial dimensionality, d. For d = 3, we estimate Pc > 3- 
The phase diagram in the plane spanned by the disorder, 
A, and the (mean) AFM coupling constant, J, for d < 4 
and p < Pc is schematically depicted in Fig. |^. 
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FIG. 1. Schematic phase diagram for d < 4 and p < pc, 
showing the antiferromagnetic (AFM) and random-singlet 
(RS) phases in the J — A plane. For A = 0, one has an 
AFM for J > Jc, and a Fermi liquid for J < Jc- The depen- 
dence of the phase boundary on the fluctuation parameter u 
is not shown. See the text for further explanations. 
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In agreement with the conclusion of Ref. one is 
in the RS phase for sufficiently small values of J for all 
A > 0, and also for sufficently small A > for all values 
of J. However, for Ai < A < A2 there is AFM LRO for 
sufficiently large J. The transition to this AFM state is 
non-Gaussian in nature, i.e. it has non-mean field expo- 
nents. Still larger disorder destroys this AFM phase 0. 
At A = 0, there is the Gaussian transition described by 
Hertz from a Fermi liquid ai J < Jc to an AFM at 
J > Jc- This AFM phase is unstable against arbitrarily 
small amounts of disorder. 

For 4 < d < d* , where d* depends on p, and p < Pc, 
the phase diagram still looks like the one shown in Fig. ||, 
except that the Gaussian AFM phase extends from the 
J-axis to finite values of A, up to a value A* < Ai. For 
d > d* the disorder window Ai < A < A2 disappears, 
and the transition is always Gaussian. The non-Gaussian 
transition also disappears for p > Pc- If d > 4, one then 
has only the Gaussian transition at small A, while for 
d < 4 there is no LRO phase. 

The qualitative features of this phase diagram are due 
to the nature of the mechanism that restores the non- 
Gaussian transition, viz. quantum fluctuations that can 
lead to long-range spin correlations which quench the 
LM. In d > 4, disorder is much less destructive for LRO 
than in d < 4. For small disorder, this allows for a Gaus- 
sian fixed point where both the quantum fluctuations and 
the disorder flow to zero. For larger disorder, quantum 
fluctuations are necessary to restore LRO even in d > 4, 
and the fixed point is non-Gaussian. In d < 4, the lat- 
ter mechanism is the only one, so there is only a non- 
Gaussian transition. The qualitative dependence of the 
phase diagram on the number of spin components follows 
along the same lines: With increasing p, as with increas- 
ing d, the AFM phase shrinks. For p > pc in d < A, LRO 
can no longer be restored and one has a RS phase for all 
values of J. In d > 4 for p > pc there still is a Gaus- 
sian transition, while the non-Gaussian one disappears if 
either p or d become too large. 

We have performed a perturbative renormalization 
group (RG) calculation that corroborates the above ideas 
and results, and also sheds additional light on the struc- 
ture of the phase diagram. We also have obtained quan- 
titative results in the framework of a double e-expansion, 
working in 4 — e space and imaginary time dimensions 
1^ . To one- loop order we obtain for the correlation length 
exponent 1/ and the dynamical critical exponent z at the 
non-Gaussian transition, 
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— [3pe + 4(p + 2)e.] , (la) 



^-2+^^— ^[(4-p)e + 4(p-t-2)e.] . (lb) 

To this order, the exponent rj = 0. All other static expo- 
nents can be obtained from the usual scaling laws pO|. 



In what follows we describe these explicit calculations, 
and show how they lead to the above conclusions. 

Our starting point is Hertz's action for an itiner- 
ant quantum antiferromagnet, which is a (/)^-theory for 
a p-component order parameter field (f> whose expecta- 
tion value is proportional to the staggered magnetization. 
The bare two-point vertex function reads. 



r(q, uj„ 



t + cf 



(2) 



Here t denotes the distance from the critical point, q is 
the wavevector, and w„ is a Matsubara frequency. We 
modify this action by adding disorder in the form of a 
'random mass' term, i.e. we consider t a random function 
of position with a Gaussian distribution with mean t and 
variance A. We use the replica trick [|ll[ to integrate out 
this quenched disorder and obtain an action, 

d^dy f^'^ drdr' ^ 0"(x,r) 



xr(x-y,r-T') r(y,r'^ 
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Here a and [3 are replica indices, r denotes imaginary 
time, r(x, t) is the Fourier transform of the vertex func- 
tion given in Eq. (|2), and u is the coupling constant of 
the usual ^''-term Q. Note that m, A > 0, so the pres- 
ence of disorder (i.e. A ^ 0) has a destabilizing effect on 
the field theory. 

Let us first reconsider the clean case, A = 0. We define 
the scale dimension of a length i to be [L] = —1, and 
that of time to be [r] = —z, with z the dynamical critical 
exponent, and look for a Gaussian fixed point where z — 
2 and 77 = 0. Power counting in d dimensions shows 
that the scale dimension of u is [u] — A — (d + z), so u is 
irrelevant for all d > 2, and the Gaussian FP is stable 
In contrast, the term ^ A carries an extra time integral, 
so with respect to the Gaussian FP we have [A] = 4 — d. 
Hence the disorder is relevant for d < 4, and the Gaussian 
FP is no longer stable in the presence of disorder. This 
instability of the Gaussian FP can also be inferred from 
the Harris criterion p^ . 

In order to see whether there is any other FP that 
might be stable instead of the Gaussian one, we have per- 
formed a one-loop RG calculation for the model, Eq. (|^). 
A simple momentum-frequency shell calculation yields 
the following flow equations. 
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Here I = \\ih with b the RG length scale factor, and we 
have scaled u — > u/24, and A — ^ A/2. Our model is 
formally very similar to the classical model studied in 
Refs. Q, and our flow equations, Eqs. (||), can easily be 
mapped onto theirs. A controlled loop expansion requires 
a double expansion in two small parameters, which in the 
present case take the form of e = d — 4, and the number 
of time dimensions e,-. The physical case is, of course, 
e — — ^1 a-iid the expansion in is probably ill- 
behaved 1^ . It is therefore reassuring to see that the FP 
structure of the flow equations does not change if one 
formally sets e-r = !> see the discussion below. 

Equations (jj) allow for four FP, which we denote by 
(M*,A*,i*). To determine the stability of the critical 
surface, it suffices to discuss Eqs. (Hi©. The four 
FP are: (1) A Gaussian FP with m* = A* = 0, (2) 
a Random FP with u* = 3(e + 4e^)/2(p - 1), A* = 
[(4 - p)e + 4(p + 2)e^]/32(p - 1), (3) an Unphysical FP 
with u* = 0, A* = -e/16, and (4) a Classical FP with 
u* = 6(e - 2e^)/(p + 8), A* = 0. The two FP that 
correspond to critical points in the present context are 
the Random and the Gaussian one. A linear stability 
analysis within our one-loop calculation shows that the 
Gaussian FP is stable for d > 4, and unstable for d < 4, 
independent of p. The Random FP is stable provided 
that 3pe > 4(p — 4)eT-. It thus is always stable for d < 4 
and p < 4, which includes the physical case c? = 3, p = 3. 
For e — €r, large values of p will destroy the stability of 
the Random FP in d < 4, while for sufficiently small p it 
is stable not only for d < 4, but also in a range of d > 4. 
For d > 4, the Random and the Gaussian FP thus may 
both be stable. The eigenvalues for the Random FP are 
complex, so the corrections to scaling at this transition 
are oscillatory in nature The remaining two FP do 
not describe phase transitions. For d < 4, the Unphys- 
ical FP is inaccessable from physical, i.e. positive, bare 
values of A. For d > 4 it is unstable. The Classical 
FP is the usual stable (for d < 4) FP for a clean, clas- 
sical (cr = 0) system. For a clean quantum system, it 
becomes unstable since the dynamical exponent z low- 
ers the upper critical dimension, and the Gaussian FP is 
stable instead. The Gaussian FP then in turn is unsta- 
ble against disorder. By linearizing Eq. (|^) about the 
Random FP, we obtain the correlation length exponent, 
Eq. (|l^). Eq. ( p^ ) and 77 = follow from the renormal- 
ization of the terms uin and q^, respectively, in Eq. (|^): 
The latter is not renormalized to one-loop order, and the 
renormalization of the former yields z = 2-\- 4A*. 

In order to determine the regions of attraction for the 
FP, we have solved the flow equations numerically for var- 
ious values of e and e^. The results did not qualitatively 



depend on the precise values as long as 2eT > e. We 
restrict our discussion to the physical quadrant w > 0, 
A > 0. For d < 4 and p > pc neither FP is stable, 
and one finds runaway flow for all initial conditions with 
A ^ 0. This runaway flow we interpret as indicative of 
LM formation, with the scale be ^ A^^/^ at which A di- 
verges the (mean) extension of a LM. We have obtained 
further support for this interpretation by constructing a 
classical instanton solution for the unstable fleld theory 
[Q. The physical meaning of the instanton is a LM, and 
be sets the scale in the instanton equation. Application 
of the methods of Refs. Q and Q to these LM will then 
yield a RS phase. 




FIG. 2. Schematic flow diagram on the critical surface for 
p < Pc and d > 4 (but d < d* , see the text). The dotted line 
separates the region of attraction for the Random FP (R) 
from a region of runaway flow. Both the Random FP and the 
Gaussian FP (G) are stable, while the Unphysical FP (U) is 
unstable. The broken line is a separatrix that marks the limit 
of the region of attraction of the Gaussian FP. 

For p < Pc (with pc = 16eT-/(4er — 3e) in one-loop 
approximation), the Random FP is stable and has a fi- 
nite region of attraction. With increasing distance of 
the initial values from the FP we still found the flow to 
be attracted by the FP, but only after oscillatory excur- 
sions to large values of A and u, where the one-loop ap- 
proximation is no longer valid. It is natural to interpret 
these large excursion as runaway flow indicative of the RS 
phase. This interpretation gets support from the two- 
loop flow equations derived by Boyanovsky and Cardy 
Q. Solving these equations for d < 4, with a reinterpre- 
tation of the parameters as appropriate for the present 
problem, we have found that the region of attraction for 
the Random FP is flnite, with a limit cycle separating 
the region of attraction from a region of runaway flow. 
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We conclude that the flow diagram on the critical surface 
(t = t*) is quahtatively as shown in Figs. ||, |3[ For d > 4, 
both the Gaussian FP and the Random FP are stable. 
The former is attractive for small values of A, and its 
region of attraction is limited by a separatrix that ends 
in the unstable Unphysical FP, see Fig. ||. For disorder 
values above that separatrix the flow is qualitatively the 
same as in d < 4. These flow diagrams correspond to 
the phase diagrams discussed above and shown for c? < 4 
in Fig. 1^. We finally discuss the general structure of 
the phase diagram, in particular why LRO occurs only 
in a disorder window. As we have seen, the nontrivial 
structure of the phase diagram arises from competition 
between fluctuation induced LRO, and LM. Technically, 
the instanton or LM contribution to the free energy goes 
like — A", with a some positive exponent. For small bare 
values of A, the instanton scale be '-^ A~^/' is large, and it 
takes many RG iterations to reach it. The renormalized 
A is therefore large, and the free energy gain due to the 
formation of LM is larger than that due to forming LRO. 
The same is true for large bare values of A. However, for 
intermediate disorder values LRO is energetically favor- 
able. The approach to the AFM FP then leads to a large 
correlation length, i.e. long-ranged spin-spin correlations 
that fall off only as l/r''~^+''. Such long-ranged correla- 
tions quench the LM |15|, which is why the rare regions 
discussed in Ref. 0] cannot preclude the AFM transition. 








G 







u 



FIG. 3. Same as Fig. |, but for d < 4. 



In summary, we have found that long-range order in 
3 — d disordered AFM systems can be restored by quan- 
tum fluctuations which destroy the random-singlet phase. 
Our most striking prediction is the possibility of re-entry 
into an AFM state with increasing disorder. Any at- 
tempts to check this prediction experimentally, however, 
should keep in mind that changing the disorder usually 



also changes other parameters, e.g. J and u. One might 
also try to confirm our result by looking for a transi- 
tion between an antiferromagnetic and a random-singlet 
phase with exponents that are neither mean-field like (as 
at T = in the clean case), nor classical (as at T > 
in the disordered case) ||l6| . Although the present theory 
has been formulated for itinerant electrons, this could be 
done in either metallic or insulating systems. 
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